Fluctuation Effects on Quadratic Autocatalysis Fronts 
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A Markov chain model for spatially distributed autocatalytic systems with a quadratic reaction 
rate is considered. An approximate solution for the local probability distribution is obtained in 
the form of a perturbation expansion for the regimes where diffusion is relatively fast. Using 
this approximate distribution, properties of the chemical wave fronts found in these autocatalytic 
systems are studied, and deviations of the minimum propagation velocity and the concentration 
profile from deterministic predictions are analyzed. A comparison with numerical results from 
lattice-gas automaton simulations is also provided. 
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I. INTRODUCTION 

Chemical wave fronts are simple examples of spa- 
tiotemporal patterns supported by noncquilibrium chem- 
ical systems .Em The dynamics of such fronts is gener- 
ally described by a set of deterministic reaction-diffusion 
equations; a wave front is defined as a special solution of 
these equations, unformly translating in space with con- 
stant velocity and connecting two distinct homogeneous 
stationary states. In particular, the wave front solutions 
of the reaction-diffusion equation, 



dc 

at 



V 2 c + c (1-c) 



(1) 



have been studied extensively. Here c is the (scalar) or- 
der parameter field (c varies between and 1), and t is 



the time variable. Equation 
FisherO and Kolmogorov et al. 
ulation dynamics and has since found applications in a 



was first considered by 
in the context of pop- 



variety of .other disciplines, including nonlinear chemical 
kinetics.tro From a chemical viewpoint, equations of 
this type describe the evolution of a spatially distributed 
autocatalytic system with a quadratic reaction rate (re- 
ferred to in this paper as a quadratic autocatalysis sys- 
tem). 

It is well known that, given appropriate initial condi- 
tions, eq. (Q) admits (one-dimensional) wave front so- 
lutions of the form c = c(x — vt), where v is the front 
velocity and x is the spatial coordinate in the direction 
of front movement. The front velocity must not be less 
than the minimum-, welocity v mm = 2 in ordepdor the 
wave front to exist.tm Further analysis (see refill for de- 
tails) shows that the solution with v = w m in is marginally 
stable and, therefore, any wave front initiated with a suf- 
ficiently steep initial profile will eventually relax to this 
minimum-velocity solution. No exact analytic solution 
of eq. (Q) is known for general v; however, it is possible 
to derive the following approximate form fqs the front 
profile, which is uniformly valid for v > v m infj 



c(z) = 



1 

1 + e z / v 



P z /v 



ez/vY 



In 



Ae z / V 



(1 + e z / v )' 



(2) 



where z — x — vt is the spatial coordinate of the reference 
frame moving with the wave front. 

This deterministic description, although capable of 
capturing the gross features of the front dynamics, is 
less successful with respect to finer details. For instance, 
reaction-diffusion equations do not provide an adequate 
framework for understanding the influence of fluctuations 
intrinsic to the local reactive and diffusive dynamics on 
the macroscopic properties of the wave front. This is due 
to the fact that reaction-diffusion equations are essen- 
tially mean-field equations in which all correlations built 
by the local fluctuating dynamics are neglected. How- 
ever, recent numerical studies on chemical fronts have 
shown that fluctuations significantly affect both the prop- 
agation velocity and the concentration profile of the wave 



front 



In this paper we present a theory which systematically 
accounts for the effects of fluctuations on the properties 
of the wave fronts in quadratic autocatalysis systems. 
Our analysis is based on the perturbation theoretic for- 
malism developed earliec-and applied to an oscillatory 
reaction-diffusion system. I 1 5 The system is described by a 
Markov chain modelE3E3 whose mean-field kinetic equa- 
tion is the generalized Fisher-Kolmogorov equatiorL_Wc 
use a Chapman-Enskog-like perturbative techniques to 
extract the corrections to the mean-field local probabil- 
ity distribution due to fluctuations in the local dynam- 
ics. The formal expansion parameter 7 gauges the ratio 
of characteristic time scales of diffusion and reaction pro- 
cesses, and is small when diffusion is fast enough, i.e. the 
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deviation from the mean-field behaviour is small. The ki- 
netic equation for the local concentration obtained from 
the generalized local distribution contains a reaction rate 
term which is not quadratic in the concentration. We 
further analyze the wave front solutions of the general- 
ized kinetic equation and determine the dependence of 
the front properties on the diffusion coefficient. We also 
compare these results with the deterministic predictions 
and the data obtained in lattice-gas automaton simula- 
tions of the quadratic autocatalysis system. 

The paper is organized in the following manner. 
In Sec. O we introduce the Markov chain model for 
quadratic autocatalysis systems and review the mean- 
field results for the properties of the wave fronts in this 
model. In Sec. JTTT]- using the perturbative technique de- 
veloped in Ref . ,E3the wave fronts are analyzed in a frame- 
work that goes beyond mean-field theory. The depen- 
dence of the minimum propagation velocity and the front 
shape on the diffusion coefficient, including the effects of 
fluctuations, is derived and compared to the mean-field 
results. Section IV presents the results of lattice-gas au- 
tomaton simulations of a quadratic autocatalysis system 
and compares these results with the predictions of the 
perturbation theory. Finally, Sec. [y| contains a discus- 
sion of the results of our study. 



II. MARKOV CHAIN MODEL 

We consider a system which consists of a chemical 
species, A, diffusing and reacting in solution. We assume 
a discrete-space, discrete-time description of the system 
in which space is partitioned into cells of fixed volume and 
time into intervals of unit length. The cells can be alter- 
natively viewed as nodes of a lattice; reactant particles 
residing at a node may react according to the reaction 
mechanism, with probability determined by combinato- 
rial rules, or execute diffusive jumps to the neighboring 
nodes. Furthermore, we impose exclusion principle, ac- 
cording to which a node cannot be occupied by more than 
N reactant particles. 

The following "particle-hole" reaction mechanism is 
used in this model: 

A + A* h 2A, 
Ah%. 

Here k\ and k 2 are rate coefficients and an asterisk de- 
notes a vacancy, i.e. an empty space at a node that can 
be occupied by a particle of A. The concentrations of 
pool chemical species (i.e. species whose concentrations 
are fixed by constraints) are incorporated in the rate con- 
stants. Reactions of this type that depend on the con- 
centrations of both particles and particle vacancies are 
common in surface chemistry and in biochemistry. The 
solvent particles are assumed to be chemically inert. At 
any time the state of the system is completely specified 
by the number of the particles of A occupying every node. 



In the following, we set k\ = 1; one can easily show that 
this is always possible by appropriate rescaling of time. 
The probabilities per unit time of the two reactions, de- 
termined by the combinatorics, are, respectively, 



Pi = 



N -1 



a(r) (N - a(r)) , p 2 = hk 2 a(r) 



Here ft. is a parameter that sets the time scale of the 
Markov chain, and a(r) is the total number of particles 
at the node with position vector r. The factor (N — 
is included so that the mean-field kinetic equation (i.e. 
reaction-diffusion equation) for the concentration of A 
can be written in a neat analytic form (see discussion 
below) . 

The inert solvent particles serve to randomize of the 
diffusive jumps that the reactant particles execute. The 
probability that a particle at the node with position vec- 
tor r will execute a diffusive jump within a unit time in- 
terval is equal to the ratio a(r)/N. The direction of the 
jumps is selected at random for each time interval and 
is the same for all nodes in the system. One can show 
that this diffusion mechanism yields the correct form of 
the diffusion equation for the local concentration. 

The Markov chain which describes evolution of the 
full probability distribution function is constructed from 
the successive application of the reaction and diffusion 
Markov chains. The numerical simulations reported in 
Sec. ^ utilize this full Markov chain, without any fur- 
ther assumption or simplification. 



A. Local Markov chain dynamics 



To facilitate theoretical analysis of the Markov chain 
model an additional assumption is employed. Numerical 
studies of the diffusion Markov chair£§E3 showed that for 
systems with spatial dimensions two or greater the dif- 
fusion rule described above leads to rapid de-correlation 
of the particle number densities at different spatial lo- 
cations. Consequently, the full probability distribution 
for the entire system may be factored into a product 
of single- node probability distributions. In the follow- 
ing, we assume that the system's dimensionality is large 
enough so that this factorization approximation is valid 
and, for sufficiently long times, the de-correlation allows 
description of the evolution in terms of the local proba- 
bility distribution function P(a(r),n). 

The evolution is due to two competing processes, 
namely local reactive events and diffusive transport of 
particles. Correspondingly, the evolution equation for 
the local distribution function P(a(r),n) can be written 
in the following form: 
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P(o(r), n + 1) - P(a(r), n) = (>yW R + W D ^j P(a(r), n). 



(3) 



Here n is the (integer) time variable and W D and W R 
are the evolution operators for pure diffusion and reaction 
processes, respectively. The dimensionless parameter 7 
gauges the relative contributions of these two processes 
to the overall dynamics. A natural choice for the defini- 



tion of 7 is the ratio of characteristic time scales of the 
diffusion and reaction processes. 

The matijix elements of the diffusion evolution opera- 
tor W D areB 



W a{r),a'(v) 



X(r,n) ( o'(r) 



N 

x(r, n. 
N 



1 



N 

a'(r) 
N 



I d /( r ), a(r )_i + —jy- II - 
a '( r ) I X( r ,n) 



N 



N 



X(r,») 
N 

<V(r),a(r) 



<5a'(r),a(r) + l 



(4) 



Here x( r ; n ) is the mean particle number density (i.e. 
concentration) averaged over the immediate neighbour- 
hood of node r, i.e. 



X(r,n) = — V a(r', 

711 ' * 



(5) 



r'eTV(r) 



where m is the coordination number of the lattice, AA(r) 



is the neighbourhood of node r and a(r', n) is concentra- 
tion at node r' at time n, 

a(r',n)= ^>(r')P(a(r'),n) . 

a(r') 

One can show that the stationary local distribution for 
the W D operator alone is binomial, 



N \ fa(r) 



a(r)J V N 



o(r) 



o(r) 

N 



jV-a(r) 



(6) 



The matrix elements of the reaction evolution opera- tor W R are 



W R , s lM = 



N ■ 



- a'(r)(N - a'(rfj S a , {r)Mr) _ 1 + hk 2 a'(r) 6 a , [rha{r)+1 

— a'(r)(N-a'(r)J + hk 2 a'(r) 



(7) 



Note that the operator W satisfies the exclusion prin- 
ciple by construction since the probability of creating a 
particle of A at a fully occupied node vanishes. 

If &2 > 0, the stationary distribution for the pure re- 
action Markov chain is 8 a r T ) o. This strongly correlated 
distribution is the consequence of the absorbing bound- 
ary at a(r) = inherent in the pure reaction dynamics. 
Indeed, as one can see from (Q), once a system evolving 
under the pure reaction Markov chain reaches the state 
with a(r) = (i.e. totally empty), it remains trapped 
in that state for all future times. In the special case 
k% = Q, there is an additional absorbing boundary at 



a(r) = N, hence the stationary distribution has the form 
n< ^a(r),o + ^a(r),Afj where n and / are constants satisfying 
n + I — 1 and determined by the initial conditions. 

In a spatially-distributed system with finite diffusion 
coefficient, the analytic form of the local probability dis- 
tribution is determined by the interplay between the re- 
active and diffusive processes described by the operators 
W and W D . Namely, as the diffusion coefficient in- 
creases, the local distribution tends to the uncorrelated 
binomial form prescribed by the pure diffusion process, 
with time dependence incorporated trivially through the 
mean particle density, viz 
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P B (a(r),n) 



N 
a(r) 



a(r, n) 
N 



a(r) 



a(r, n) 
N 



N-a(r) 



(8) 



Conversely, as the diffusion coefficient tends to zero, cor- 
relations characteristic of the stationary distribution of 
the pure reaction Markov chain contaminate the bino- 
mial form (|J) to an increasing degree. The subject of 
this paper is the effect of these correlations on the prop- 
erties of the macroscopic structures, such as wave fronts, 
existing in the system. 



B. Mean- field dynamics 

Consider the mean-field dynamics of the reaction- 
diffusion Markov chain. The mean-field local distribution 
function is the time-dependent binomial (^); the kinetic 
equation for concentration of A can be easily obtained 
using it and eq. (0); after some analysis, we find 



N 

i(r, n + 1) - a(r, n) = ^ a(r) (w D + jW R ^j P B (a(r),n) 



a=0 
1 



Aa(r, n) + 7/1 a(r, n) /3 — 



mN 

D Aa(r, n) + jR(a(r, n)). 



a(r, n) 
N 



(9) 



Here j3 = l — k2,D = (mN)^ 1 , A is the discrete Lapla- 
cian operator and R(a(r, n)) denotes the mass-action law 
terms. Introducing the rescaled variables, 



h — > 7/3/i, a(r, n) 



a(r, n) 

N/3 



and the length unit, L = \J h/D, one can rewrite eq. (j^) 
as follows: 



a(r, n + 1) — a(r, n) = Aa(r, n) + h a(r, n) ^1 — a(r, 71)^ . 



(10) 



Clearly, eq. ( fio| ) is just a discretized form of eq. (0). 



III. PERTURBATION THEORY FOR THE 
WAVE FRONT SOLUTIONS 

We now turn our attention to the analysis of the evo- 
lution equation (^). We derive an approximate ki- 
netic equation for the local concentration field valid 
for small 7 and analyze the properties of its wave 
front solutions. 



A. Local concentration dynamics for 7^1 

We seek the solution of eq. (||) in the form of a regu- 
lar perturbation series with time-dependent binomial 
distribution as the leading term, 



P(a(r),n) = P B (a(r),n) + V 7 fc P fc (a(r),n). (11) 



fe=i 



In the following, we will keep only the leading and first- 
order terms of the series ([ll ) . The starting point of our 
analysis is the evolution equation for the first-order term, 
Pi(a(r),n\ which is obtained by substituting eq. fill ) 
into eq. (k|), collecting terms of similar order in 7 and. 
using the solvability conditions (cf. eq. (20) of Ref.tij 
and discussion therein), viz 



Pi(a(r),n + 1) -P x (o(r),n) = W D P 1 (a(v), n) + SP B (a(r), n), 



(12) 
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where the operator S is defined as follows: P 1 (a(r) n) = — ( W D \ * SPB(a(r) n) (14) 



S = W R -R(a(r,n))-^ . (13) x 

^ a where I VF 15 J is the operator reciprocal to W^ . Using 

Suppose that P\{a{v),n) changes slowly with time. (|l4|) and (pj|), one can now determine the dynamics of the 

Then, we can write the formal solution of eq. (11^) as local concentration and we obtain 



3(r,n + l) -a(r,n) = ^ a(r) Uw R + W D ^ (p B (a(r),ri) + -/P 1 (a(r),' l 

o(r)=0 

1 N -1 

— Aa(r,n)+ 7 i?(a(r))- 7 2 £ a(r) (#°) SP s (a(r),n). (15) 



a(r)=0 



( |l5| ) is presented in the Appendix. We find 

The calculation of the 0(7 2 ) correction term in eq. 



J2 a(r)^(^) _1 ^ B (a(r),n) = ^a(r,n)(l-^) 2 . (16) 

a(r)=0 ^ ' 

The final form of the kinetic equation for a(r, n) is ob- over, at this point we can set the ordering parameter 7 



tained by substituting (|16|) into (15) and transforming to unity. In this way we come to 



the concentration variable, a(r,n) — ► a(r,n)/N. More- 



/ \ Nh 1 \ z 

a(r, n + 1) — a(r, n) = D Aa(r, n) + h a(r , n) \J3 — a(r , n) J — — — - a(r,n) yl — a(r,ri)j , (17) 

where = 1 — fe. Note that eq. ( fl7| ) can be derived from lows). 

the exact solution of ( |l2| ) by an asymptotic expansion for For any fixed h, eq. (^) can be viewed as a discretized 
small 7 of the integral terms in-the corresponding kinetic version of the following partial differential equation (writ- 
equation (cf. eq. (25) of Ref.tS and discussion that fol- ten here for the case of a planar front): 



da(x,t) jjd 2 a(x,t) | _ 



dt dx 2 



+ a{x,t)(B-a{x,t^-Ka{x,t)(l-a(x,t)\ , (18) 



where D = (mNh)" 1 , K = -f^. Note that, as defined, Equation (jT|) can be written in terms of the same vari 



K is inversely proportional to the diffusion coefficient D, ables as the mean-field reaction-diffusion equation (10) 
K = a/D, where a = — ■ with a rescaling of space, x — > x/y/D, as follows: 



da(x,t) d 2 a(x,t) 
dt = dx 2 



+ a{x,t) (p-a(x,t)j - a(x,t) (l-a(x,tf) . (19) 
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The reaction rate in eq. ( |l9| ) is not quadratic in a(x,t), 
in contrast to eq. (|l0|), because of the 0(D~ 1 ) correc- 
tion term. For convenience, all further analysis will be 
carried out using the continuous-space, continuous-time 



equation (19) rather than the finite-difference equation 
(|l7|). Since eq. (17) is obtained from eq. ( |l9|) by Euler 



discretization and a linear transformation of variables, 
both equations describe essentially the same dynamics. 



B. Analysis of the wave front solutions 



Following the standard formalism,!!!' i we recast eq. 
( |l9| ) in terms of the variable z, the coordinate of a refer- 
ence frame moving with the wave front, i.e. z = x — vt, 



d 2 a 



v— = a (a — (3) 
az 



a — — \2 

— a (1 - a) 



(20) 



Here we will consider only wave fronts moving to the right 
(v > 0); generalization of our results to the left-moving 
fronts is trivial. Letting p = da/dz, we can rewrite eq. 
(EG) as a system of first-order ODEs, 



zone II: the point at the origin is unstable, the point 
at Tiis is unstable and moves away from the 
origin as D increases (ai s > 0), the point at 
a,2 S is stable and moves toward the origin as D 
increases (a2 S > 0), 



P- 



(21) 



-vp - 



■0) 



all 



where the dot indicates differentiation with respect to 
z. This system posesses up to three fixed points, one 
at the origin of the phase plane and the other two at 
(a,p) = (ai s ,0) and (a2 S ,0), where 



D - y/D 2 +4a(/3- 1)D 

ii , — 1 — , 

2a 

( _ D + y/ D 2 + 4a(/3 - 1)D 
2a 



a.2 S = 1 - 



(22) 



zone III: the point at the origin is unstable, the point 
at a\s is stable and moves toward the origin as 
D increases (au < 0), the point at a,2 S is un- 
stable and moves away from the origin as D 
increases (a2s < 0), and 



zone IV: the point at the origin is stable, the point 
at ai s is unstable and moves away from the ori- 
gin as D increases (a± s > 0), the point at a,2s 
is unstable and moves away from the origin as 
D increases (a2s < 0). 



From a linear stability analysis, we find the follow- 
ing relation for the eigenvalues associated with the fixed 
points: 



a 
D 



(1 -a) (1 - 33) - 2a + = 0, (23) 



where a = 0,ai s , or a,2s- The diagram in Fig. |l| sum- 
marizes the behaviour of the fixed-point structure with 
respect to D and f3, as determined by eqns. ( p2| ) and 
(p3[). The four zones in the diagram correspond to dif- 
ferent linear stability properties of the fixed points as 
follows: 

zone I: only the fixed point at the origin exists and is 
unstable, 



The two lines that define these zones are: a) D = 
4a(l — (3), the line of parameter values at which the fixed 
points at au and 52s first emerge as a stable, degenerate, 
zero-eigenvalue node (dashed line), and b) D = a/ (3, the 
line where the fixed point at the origin changes its sta- 
bility (solid line). We observe that formation of a wave 
front replacing the steady state at the origin with the 
steady state at 3i s is possible if the diffusion coefficient 
is greater than the critical value D cr — a/ [3. Indeed, in 
that region of the diagram the fixed point at the origin 
is stable, the point at 3i s is unstable (with a\ s > 0), and 
a2s is negative. Moreover, it can be easily seen that 3i s 
tends to (3 as D increases. Thus, in the limit D — > oo 
the fixed-point structure of (21) is equivalent to that of 
eq. © of Sec. H save for additional fixed point in an 
unphysical area of the phase plane. 



G 



CO. 




0.02 0.06 
D 



FIG. 1. Diagram of the fixed-point structure of the ODE 
system (21) as a function of D and /3, with a = . The crit- 
ical diffusion coefficient line D CI = a/p is indicated by a solid 
line. The dashed line denotes the locus of (D, (S) values where 
the movable fixed points first emerge as a degenerate node. 
For description of the fixed-point behaviour in the different 
zones in the diagram, see text above. 

In order to represent a physically meaningful state of 
the system, a wave front solution of (21) with < (3 < 1 
must exhibit no oscillatory behaviour near the leading 
edg e of the front. This implies that both roots of eq. 
( p3| ) with a = should be real. One can show by elemen- 
tary analysis that this requirement is satisfied if 



before the measurements were taken to ensure that the 
wave front relaxed to its stable, uniformly propagating 
form. In our simulations, the period of transient dynam- 
ics was found to be relatively long (~ 10 5 automaton time 
steps). Within this period, the front width was observed 
to grow diffusively; however, no significant growth in the 
front width was detected after the transient regime. 

In Fig. H we compare the front velocity observed in 
these simulations as a function of the diffusion coefficient 
D with the minimum propagation velocity predicted on 
the basis of the present theory (eq. (Q)) and the mean- 
field, reaction-diffusion equation (I1Q). 




(24) 



For any finite D the minimum propagation velocity given 
by ( |24| ) and the concentration at the upper plateau of the 
wave front (ai s in eq. (22)) are both less than their re- 
spective values predicted by the mean-field theory and 
converge to those values as D — > oo. 



IV. NUMERICAL RESULTS 



FIG. 2. Dimensionless front velocity v versus the diffusion 
coefficient D. Circles denote the data obtained in lattice-gas 
automaton simulations. The minimum propagation velocity 
given by the present theory (eq. (p^)) is plotted with a solid 
line. The results of the mean- field theory are indicated by a 
dashed line. Parameters a and /3 are equal to and 0.667, 
respectively. 

We observe that the propagation velocity predicted by 
the present theory is in good agreement with the data 
from the numerical simulations. 



In order to test the accuracy of the results obtained 
in the preceding section, we numerically solved the full 
Markov chain dynamics using a lattice-gas automaton 
algorithm. The simulations were performed on a thin 
(4 x 3000 nodes) strip of a 2-dimensional triangular lat- 
tice; wave fronts were generated from a step-like initial 
concentration profile placed in the middle of the lattice. 
The value of the rate coefficient ki used in these sim- 
ulations was 0.333, so that (3 — 0.667. The exclusion 
parameter was fixed at N — 21 which, together with 
the coordination number for triangular lattice (m = 6), 
gives a = j^. The value of the diffusion coefficient 
D = (mNhy 1 varied between 0.1 and 1.0. All of the 
front properties were measured using the mean concen- 
tration profile obtained by averaging over 10 realizations 
of the dynamics. Sufficient time was allowed to elapse 



V. CONCLUSIONS 

The generalized kinetic equation derived in this study 
admits travelling wave solutions whose properties are dif- 
ferent from the predictions of mean-field theory. There 
is a non-zero critical value of the diffusion coefficient be- 
low which no stable wave front solution exists. More- 
over, the concentration at the upper plateau of the wave 
front varies with the diffusion coefficient as well as the 
kinetic parameters of the system. Most significantly, the 
(dimensionless) minimum front velocity decreases as the 
diffusion coefficient is decreased. These effects have their 
origin in the local correlations built by the fluctuating 
reactive dynamics which persist for finite values of the 
diffusion coefficient. 
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For finite values of the diffusion coefficient, measure- 
ments of the front velocity in two-dimesional lattice-gas 
simulations of the full Markov chain model confirmed 
that the present theory describes the departure from the 
mean-field results with quantitative accuracy. Similar 
simulations were also carried out on one-dimensional sys- 
tems. The results of these simulations showed the same 
qualitative trend as in two dimensions. However, because 
the primary assumption of our theory (the factorization 
of the full probability distribution) is not valid for one 
spatial dimension, the quantitative agreement between 
theoretical data and the simulation results is much poorer 
than that for two dimensions. 



There is a number of numerical studiest3iiSI^I on the 
breakdown of mean-field descriptions for quadratic au- 
tocatalysis fronts using a variety of microscopic models, 
none of which is the same as that employed in our study. 
The. lattice-gas automaton simulations of Lemarchand et 
al$3 have shown that the mean-field results remain valid 
for the range of diffusion coefficient values used in their 
study. Later investigations using a Langevin equation ap- 
proach reported an increase in the front velocity above 
the mean-field value.HEj Monte Carlo studies of Rior- 
dan et al. have focused on the front width and have 
found that the width grows algebraically in one and two 



space dimensions, while in higher dimensions the width, 
dynamics is adequately described by mean-field theory SB 
In our simulations on two-dimensional systems with small 
transverse size, the front width grows diffusively during 
a very long transient period but then saturates and the 
front evolves without change in width. The effects of 
microscopic fluctuations on quadratic autocatalyis fronts 
clearly present a number of interesting features. 
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APPENDIX 



To compute the 0(7 2 ) correction term in eq. (15), 
note the following properties of operator W D (the proof 
can be obtained bjf straightforward calculation, as in Ap- 
pendix A of Refill): 



N 

a(r)(w D Y 'X(a(r),n) = 

a(r)=0 
N 

a\v)(w D ) P X(a(r),n) = 

o(r)=0 



N 



N 



2_V 

' N 



a(r)X(a(r),n), 



a(r)=0 
N 



]T a 2 (r)A(a(r),n), 



(25) 



a(r)=0 



where (w D ^ is to be understood as a time-ordered 

product of operators W D taken at any p moments of 
time between and n. X(a(r), n) is any function of a(r) 
and time which has the following properties: 



N 

E 

a(r)=0 



X(a(r),n) = 0, 



N 

E 

a(r)=0 



a(r) X(a(r),n) = 0. 



(26) 



We remove one of the operators W D by adding the 
inverse operator fw^M , taken at an appropriate mo- 
ment of time, to the left of ^W^ 15 ^ in (|25|). Bearing in 
mind the property (|25|), we find 



]T a(r)^ 1 (W D ) X(a(r),n)=(--\ £ a(r)A(a(r), 



a(r)=0 
N 



a(r)=0 
„ v p-l N 

a\v)W D -' (w D ) P X(a(r),n)= (--) £ a 2 (r)A(a(r), n). 

a(r)=0 ^ ' a(r)=0 



(27) 
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Note that if (§|) holds for X(a(r),n), it also holds for 
(w°y X(a(r),n), i.e. W D is a fvmction-to-function 
map within the class of functions satisfying (p6|). Fur- 
thermore, SPg(a(r), n) belongs to that class of functions, 
as can be easily seen from the definition of the operator 
S ( |l3| ) and the following two identities: 



N _ 1 

J2 a(r) (w D ^j SP B (a(r),n) 

o(r)=0 



a 2 (r) (W D ) SP B {a{v),n) 

a(r)=0 



Using (|7|) and (p9|), after a few simple transformations, 



—l 

a{r)W R (w D ^j SP B {a{r),n) 

a(r)=0 



for the C(7 2 ) correction term. 
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